LARGE DEVIATIONS OF THE LIMITING DISTRIBUTION IN THE 
SHANKS-RENYI PRIME NUMBER RACE 



YOUNESS LAMZOURI 

Abstract. Let g > 3, 2 < r < ^(g) and ai, . . . , be distinct residue classes modulo 
, q that are relatively prime to q. Assuming the Generalized Riemann Hypothesis GRH 

' and the Linear Independence Hypothesis LI, M. Rubinstein and P. Sarnak [TT] showed 

that the vector-valued function Eq-^ai,...,a^{x) — {E{x;q,ai), . . . , E{x;q,ar)), where 
^ I E{x;q,a) = {(l){q)Tr{x;q,a) — Tr{x)), has a limiting distribution Hq:ai,...,ar which 

' is absolutely continuous on MJ" . Furthermore, they proved that for r fixed, /ig;ai,...,Or 

, tends to a multidimensional Gaussian as q ^ oo. In the present paper, we determine 

^SJ ' the exact rate of this convergence, and investigate the asymptotic behavior of the 

large deviations of ^.q-^ai,...,ar- 

h' 

^ ■ 1. Introduction 

A classical problem in analytic number theory is the so-called " Shanks and Renyi 

prime number race" (see [3]) which is described in the following way. Let g > 3 and 

^ I 2 < r < </)(g) be positive integers. For an ordered r-tuple of distinct reduced residues 

^ \ (ai, 02, . . . , Or) modulo q we denote by Pq-ai,...,ar the set of real numbers x >2 such that 

O 
p 

^ ■ Will the sets -Pg;a^(i),...,a^(^) contain arbitrarily large values, for any permutation a of 
{1,2,. ..,r}? 

A result of J. E. Littlewood [7] from 1914 shows that the answer is yes in the cases 
(g, 01,02) = (4,1,3) and (g, 01,02) = (3,1,2). Similar results to other moduli in the 
case r = 2 were subsequently derived by S. Knapowski and P. Turan [3] (under some 
hypotheses on the zeros of Dirichlet L-functions), and special cases of the prime number 
race with r > 3 were considered by J. Kaczorowski 

In 1994, M. Rubinstein and P. Sarnak [11] completely solved this problem, con- 
ditionally on the assumptions of the Generalized Riemann Hypothesis GRH and the 
Linear Independence Hypothesis LI (which is the assumption that the nonnegative 
imaginary parts of the zeros of all Dirichlet L-functions attached to primitive charac- 
ters modulo q are linearly independent over Q). To describe their results, we first define 



7r(x; g, Oi) > 7r(x; g, 02) > ■ ■ ■ > 7r(x; g, o^). 
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some notation. For any real number x > 2 we introduce the vector-valued function 

Eq-a^,...,arix) ■= {E{x; Q, tti) , . . . , E {x; Q, ttr)) , 

where 

E{x; q, a) := {(f){q)7r{x; q, a) - tt{x)) . 

Jx 



The normalization is such that, if we assume GRH, Eq.ai,...,ar{^) varies roughly bound- 
edly as x varies. Rubinstein and Sarnak showed, assuming GRH, that the function 
Eq-ai,...,ar{^) has a limiting distribution fJ.q-ai,...,ar- More precisely, they proved 

1 dx f 

(1.1) lim -/ f {Eq.au...,aAx)) — = f {Xl, ■ ■ ■ , Xr)dflq.au...,ar, 

X-^oo log X J2 X J^r 

for all bounded, continuous functions / on M^. Furthermore, assuming both GRH and 
LI, they showed that /Xg;ai,...,a^ is absolutely continuous with respect to the Lebesgue 
measure on W if r < 0(g). (When r = 0(g), fJ'q;ai,...,ar is shown to be absolutely 
continuous with respect to the Lebesgue measure on the hyperplane =0.) As 

a consequence, under GRH and LI, the logarithmic density of the set Pg-ai,...,ar defined 
by 

lim / — 

X^OologX JteP,.,a„...,arn[2,x] t 

exists and is positive. 

Here and throughout we shall use the notations ||x|| = ■\jYl'j=i x] and |x|oo = 
maxi<j<r \xi \ for the Euclidean norm and the maximum norm on W respectively. In [TT] , 
Rubinstein and Sarnak also studied the behavior of the tail /Ug;ai,...,ar(l|x|| > V^) when 
g is fixed. They showed, under GRH, that for all distinct reduced residues ai, . . . ,0^ 
modulo g we have 

(1.2) exp(-exp(c2(g)l^)) fiq-a,,...,aMM > V) <9 exp (-ci(g)v^) , 

for some Ci(g), C2(g) > which depend on g. 

In this paper we investigate large deviations of the distribution fiq-ai,...,ar uniformly 
as g — )■ 00, under the additional assumption of LI. For a non-trivial character x mod- 
ulo g, we denote by {7^} the sequence of imaginary parts of the non-trivial zeros of 
L{s,x)- Let Xo denote the principal character modulo g and define S = ^x¥=xo mod q{lx}- 
Moreover, let {f/(7^)}^^g5 be a sequence of independent random variables uniformly 
distributed on the unit circle. Rubinstein and Sarnak established, under GRH and LI, 
that the distribution ^q-ai,...,ar is the same as the probability measure corresponding to 
the random vector 



(1.3) 



^q;au...,ar = i^il, Ctl), • • • , ^(g, ^r)), 
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where 

X(,.a) = -C,(a)+ J: 

Xt^Xo 7x>0 \/i +7x 
X mod q ^ 

Xo is the principal character modulo q and 



C,{a) := -1 + Yl 1- 



b^=a mod g 
l<f)<g 



Note that for (a, g) = 1 the function Cq(a) takes only two values: Cq{a) = —1 if a is a 
non-square modulo q, and Cq{a) = Cg(l) if a is a square modulo q. Furthermore, an 
elementary argument shows that Cq{a) < d{q) 1^ for any e > 0, where d{q) = J2m\q ^ 
is the usual divisor function. 

Let CoVq-ai,...,ar t>e the covariance matrix of Xq.ai,...,ar- A straightforward computa- 
tion shows that the entries of CoVg^a^^...^^^ are 



Var(g) if j = 

Bq{aj,ak) if j ^ k, 



where 



and 



Xt^XO 7x 4 ^ 'X Xt^XO 7x>0 4 ^ 

X mod q x mod 9 

x(!)+x(f) 



Bq{a,b):= J2 E 



Xt^XO 7x>0 4 ^ 'X 

X mod q 

for (a, 6) G ^(g), where A{q) is the set of ordered pairs of distinct reduced residues 
modulo q. Assuming GRH, Rubinstein and Sarnak showed that 

(1.4) Var(g) ~ 0(g) log q and Bq{a, b) = o{(f){q) log g) as g — > oo, 

uniformly for all (a, b) G A{q). Combining these estimates with an explicit formula for 
the Fourier transform of fJ^q-ai,...,ar (s^e (3.1) below) they established, under GRH and 
LI, that 

/i,;a„...,a. (||x|| > Av^Va^) = {2ny/' [ exp f^^l_±^\ rfx + o,(l„ 

^ ^ -'l|x||>A V ^ / 



for any fixed A > 0. 

We refine their result using the approach developed in [6]. More precisely, we prove 
that the asymptotic formula (1.5) holds uniformly in the range < A < i/Ioglogg with 
an optimal error term 0^.(1/ log^ g). 
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Theorem 1. Assume GRH and LI. Fix an integer r > 2. Let q be large and 

be distinct reduced residues modulo q. Then, in the range < A < y/log log q we have 

/i.;a„...,a. (||x|| > X^^^) = (2^^^' I exp ( ^l + '-' + ^'A ^^^Qr ( -\-) 

^ ^ iiix||>A V 2 y viog qj 

Moreover, there exists an r -tuple of distinct reduced residue classes (ai, . . . , Qr) modulo 
q, such that in the range 1/A < X < 3/4 vue have 



>r 



/i,;a„...,a. (||x|| > Av^Va^) - (27r)-'/2 f exp fd±:_±^\ 

^ ' -'l|x||>A V ^ / 

Since Var(g) ~ 0(g) log g, it follows from Theorem 1 that 

(1-6) /^,;a,...,a.(||x|| > 1^) = exp [-^T^. (1 + o(l)) 

V 20(g) log g 

in the range (0(g) log g)^/^ V < (1 + o(l))(0(g) loggloglogg)^/^. Exploiting the 
results of H. L. Montgomery and A. M. Odlyzko [lO], we prove that a similar behavior 
holds in the much larger range (0(g) logg)^/^ -C <^ 0(g) logg. 

Theorem 2. Assume GRH and LL Fix an integer r > 2 and a real number A > 1. 
Let q be large. Then for all distinct reduced residues ai, . . . , modulo q, we have 

exp -C2{r,A)—— < /ig;ai,...,a,(| |x| | > 1^) < cxp -Ci{r,A)- 



0(g) log gy ' " V ' 0(?)log?. 

uniformly in the range (0(g) logg) ^''^ <^V < A0(g)logg, where C2{r, A) > ci{r,A) are 
positive numbers that depend only on r and A. 

Using an analogous approach, we prove that (1.6) does not hold when V/ (0(g) log g) 
oo as g — 7- oo, which shows that a transition occurs at \^ x 0(g) logg. 

Theorem 3. Assume GRH and LL Fix an integer r > 2 and let q be large. If 
V/{(j){q) logg) — oo and V/{(f){q) log^g) — )■ as q ^ oo, then for all distinct reduced 
residues Oi, . . . , modulo q, we have 

(-^^^^^0(^^"p {''^'^m^)) > 

and 

/i........(||x|| > V) « exp (--3(r)^^^exp (^5^^^^)) , 

where C4(r) > Cs{r), and CQ{r) > C5(r) are positive numbers which depend only on r. 

In the range V/ (0(g) log^ g) — )■ oo one can prove, using the same ideas, that there 
are positive constants Cs{r) > C7(r) such that 

(1-7) ^ ^ ^ , . ^ 




exp -exp C8(r) < /ig;ai,...,a,(| |x| | > 1^) < exp -exp cr{r 
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In particular, these bounds show that the asymptotic behavior of the tail /^g;ai,...,ar(| |x| | > 
V) changes again at x (p^q) log^ q. Assuming the RH and using the LI for the Rie- 
mann zeta function, H. L. Montgomery [9] had previously obtained a similar result for 
/ii, the limiting distribution of the error term in the prime number theorem 7r(x)— Li(a;). 
His result states that 

exp (^-cioVV exp (y2TcV^^ < Hi{\x\ > l^) < exp (^-cgW exp (^V2tcV 

for some absolute constants cio > cg > 0. A more precise estimate was subsequently 
derived by W. Monach [8], namely 

(1.8) >V)= exp (^-(e~^° + o(l))y2^exp (y2nV^^ , 

where Aq is an absolute constant defined in Theorem 4 below. 

In our case, it appears that /Xg;ai,...,ar(|x|oo > V) is more natural to study, if one 
wants to gain a better understanding of the decay rate of large deviations of ^q-ai,...,ar 
in the range V/ {4>{q) log^ q) — > oo. We achieved this using the saddle-point method. We 
also note that in contrast to our previous results, r can vary uniformly in [2, 0(g) — 1] 
as g — > oo. 

Theorem 4. Assume GRH and LI. Letq be large, and 2 <r< 0(g) — 1 be an integer. If 
^ I i^'il) log^ g) — oo as q ^ oo, then for all distinct reduced residue classes 
modulo q, the tail /ig;ai,...,ar(|x|oo > V) equals 





+ Aq- log TT, 



and 



io t Ji t 

where Io{t) = X]^o(V2)^"'/'^!^ is the modified Bess el function of order 0. 

Lastly, we should also mention that in the range V > \J (t){q) logg one may allow r 
to vary uniformly, as in Theorem 4, if one is willing to replace /i5;ai....,a^(| |x| | > V) by 
/^g;ai,...,ar(|x|oo > V) in the Statements of Theorems 2 and 3. 

2. The intermediate range: Proof of Theorems 2 and 3 

In this section we investigate the behavior of the tail /ig;ai,...,ar(| k| | > V) when 
V ^ a/0(^)^oS^ long as \^/(0(g) log^ g) — )■ as g — t- oo. Recall that Hq-ai,...,ar is 
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also the probability measure corresponding to the random vector Xq.^ai,...,ar (defined in 
(1.3)). Our idea starts with the observation that the random variables 

(2.1) YM:=xM-nxM)= y: 

Xt^XO 7x>0 V i + 7x 
X mod q ^ 

are identically distributed for all reduced residues a modulo q. Indeed, for all (a, q) — 1 
the random variables {f^(7x)}7xeS) where U{'y^) = x(a)f^(7x)) ^ire independent and 
uniformly distributed on the unit circle. Hence, if (a, q) — 1 then Y{q, a) has the same 
distribution as 

2cos(27r^(7^)) 



(2.2) y(,):= Y, E /i— . 

Xt^XO 7x>0 V4+7x 
X mod q ' 

where {^(7x)}7xeS' are independent random variables uniformly distributed on [0,1]. 
Our first lemma shows, in our range of V, that large deviations of f^q;ai,...,ar are closely 
related to those of y(g). 

Lemma 2.1. The random variable Y{q) is symmetric. Moreover, if q is sufficiently 
large, r > 2 is fixed and V > (j){q), then for all distinct reduced residues 
modulo q we have 

F{Y{q) > 2V) < l^q.,a,,...,aAM > V) < 2rP (Y{q) > . 

Proof Note that E(cxp(itcos(27r^(7^)))) = Jo(t), where Jo(t) := Em=o(-l)'"(V2)^"'M!^ 
is the Bessel function of order 0. Therefore, since the Fourier transform (characteristic 
function) of Y{q) is 

/ 



Xt^Xo 7x>0 
X mod q 



2t 



and Jo is an even function then Y[q) is symmetric. Now, 

l^q;a,,...,aMM > V) = F{\\Xq,a„...,aM > V) > P(X(g, a^) > V) = F{Y{q) > V+Cq{a,)). 
The lower bound follows from the fact that Cg{ai) < d(q) = q°^^\ On the other hand, 
if \X{q,aj)\ < Vj^ for all 1 < j < r, then ||-^g;ai,...,ar.|| < ^ ■ Hence 

l^q;a.,...,aX\\A\ > V^) < E ? (|^(g,%)| > ^) < E ^ (l^(^)l > ^ " " 

Using that V > \/0(^, |C'5(a)| = ^"^^^ for all (a, g) = 1, and Y{q) is symmetric we 
obtain from the last inequality 

l^q;a,,...,arm\ > V) < vF (^Y{q)\ > = 2rP (Y{q) > ^) , 

if q is sufficiently large. This estabhshes the lemma. 
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□ 

To investigate large deviations of Y{q), we shall appeal to the following result of H. 
L. Montgomery and A. M. Odlyzko [TO] . 



Theorem 2.2 (Theorem 2 of [lOj). Let {Yn}n>i be a sequence of independent real 
valued random variables such that E(y„) = and \Yn\ < 1. Suppose there is constant 
c > such that E{Y^) > c for all n. Put Y = J2n>i^n/rn where J2n>i < ^/ 



(2.3) P(F >V)< exp 
Moreover, if J2\r„\<T\''^'n.\~^ > 21^, then 

(2.4) P(F > 1/) > aiexp 
where ai, 02 > depend only on c. 




( 



\ 




In order to apply this result to our setting, we have to understand the asymptotic 

behavior of the sums Y^x^x^ So<7x<t V(i + t')'^' and ^x^xo ^i^>t V(i + t')- For 
a non-trivial character x modulo g, we let g* be the conductor of and x* be the 
unique primitive character modulo g* which induces x- We begin by recording some 
standard estimates which will be useful in our subsequent work. 

Lemma 2.3. Assume GRH. Let x be a non-trivial character modulo q. Then 



XI 1^ 2 = ^log9x + ^(^°slogg). 



+ 72 2 

7x>0 4 ^ 'X 



Moreover, we have 



log^x = '^(^) log? - 5Z 

; mod g \ v\l / 



and 



logp ^ , , 
2^ < log log g. 



Proof. The first estimate follows from Lemma 3.5 of p], and the second is proved in 
Proposition 3.3 of [2]. Finally, we have 



logp >r;^ logp 1 , 

V 7< V 7 + ^^ > Kloglogg, 

^P-^ ^n,2P~^ logg^ 

p\q p<(logg)^ Pk 

which follows from the trivial bound Ylp\q ^ — logg/ log 2. □ 
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From this lemma, one can deduce the more precise asymptotic 
Var(g) = (f){q) logg + O(0(g) log log g). 
Our next result gives the classical estimate for 



x/xo o<7x<r 

X mod q 



Lemma 2.4. IfT >2, then 



N,iT)= ^-^^^T log T + + O (0(g) log gT) , 

where 

" (0(g)-l)(log27r + l). 



R{q) :=0(g) (logg-E^ | 



Proof. Using Lemma 2.3 and appealing to the classical estimates for the number of 
zeros of Dirichlet L-functions (see Chapters 15 and 16 of [Ij), we get 



Xt^XO 
X mod q 

= M£l^riogT + Mr + ow,)iog,r). 

as desired. □ 
Using Lemmas 2.3 and 2.4 we establish: 

Proposition 2.5. There exists a constant To > 2 such that ifT > Tq, andlogT/ logg — )■ 
as g — )■ oo, then 

^0(g) logg log T < V V ^ < 30(g) logg log T 



X mod g 

«/g Z5 sufficiently large, and 



1 _ 0(g) logg (. p 



0(g)logg/^ , ^ /logT ^ loglogg , 1 



X mod q 

Proof. First, since logT = o(logg) then Lemmas 2.3 and 2.4 yield 

T 

^qiT) = —0(g) logg + O (0(g)Tlogr + r0(g) log log g + 0(g) logg) 

(2 5) 

= 7^0W logg 1 + O -j + — + - 

27r \ \ log g log g 1 ^ 

Hence, there exists a suitably large constant Ti > 2 such that 
(2.6) ^0(g) logg < iV,(T) < T0(g) logg. 
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if T > Ti and logT = o(logg). Now assume that T > Tf. Then 



T T - = 

X^XO 0<,,<T + Jo (i + t^) 

' mnrl /7 V V V 



X mod q 

Then, using the lower bound of (2.6) we get 

X mod q ' 

which follows from the fact that > (1/4 + for i > 1. Similarly, we obtain from 
the upper bound of (2.6) 

/■^i tN (t) tN (t) 

/ .1 ^\;,, (it«T, 0(g) log g and / dt< 0(g) log g log r. 

This implies 



J2 Yl ^ 20(g)logglogr + O(0(g)logg). 

xt^xo o<7x<r Vi + 7x 

X mod g » 



Therefore, if T > Tq for some suitably large constant Tq, then 

^0(g)logglogT< J] J] J < 30(g) log g log T. 

xt^xo o<7x<r V i + 7x 

X mod g ' 

On the other hand, we have 

^ ^ 1 _ r dN,{t) _ N,{T) , ^ r tN,{t) 

X mod q 

Thus, inserting the estimate (2.5) into the previous identity gives 



where 



XT^XO 7x>^4 ' 'X ""V4 ' - ^ " -^^ 

X mod q 



^ ^ 0(g) log g /^ log^ ^ log log g ^ ^ . 
^ T viogg logg 



Moreover, the main term of (2.7) equals 



0(g) logg f^^Qf]_ 



2nT 

as desired. □ 
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Proof of Theorem 2. Since E(cos(27r^(7^))) = and E(cos2(27r^(7^))) = 1/2 > 0, then 
we can apply Theorem 2.2 to derive upper and lower bounds for P(F(g) > V). We first 
establish the upper bound. Taking T = in (2.3) yields 

y2 \ / 



(2.8) ^{Yiq) >V)< exp , , < exp , , , , , 

^ ' \ )- i'y 32Var(g)y - 400(g) log 

if q is sufficiently large, since Var(g) ~ 0(g) log g. Now, let T = T{A) be a suitably 
large number such that logT > 40^4. Then Proposition 2.5 implies 

E E ^==^ E E ^i=>^0(g)ioggiogT>2V^, 

Xt^XO (1/4+72)1/2<2T V 4 + 7x Xt^XO 0<7x<r A / 4 + 7^ 

Xmodg ^^^q v x mod g V 

since V < A(j){q) logg. On the other hand, using Proposition 2.5 we get 

Xt^XO (l/4+72)i/2>2r 4 'X Xt^XO 7x>2r 4 'X 

Xmodg ^^^Q xmodg 

if T is suitably large. Hence, applying (2.4) we derive 

(2.9) P(F(g)>l^)>exp -Ci(A) 



0(g) log g^ 

for some positive number ci{A) which depends only on A. The theorem follows upon 
combining the bounds (2.8) and (2.9) with Lemma 2.1. □ 

Proof of Theorem 3. The result can be deduced by proceeding along the same lines as 
in the proof of Theorem 2. Indeed, using Lemma 2.1, Theorem 2.2 and Proposition 2.5, 
the lower bound in Theorem 3 can be derived by taking T = exp(5Oy^/(0(g) logg)). 
Similarly, to get the corresponding upper bound choose T = exp(V^/ (100(g) logg)). □ 



3. Approximating Hq-ai,...,ar by a multivariate Gaussian distribution: 

Proof of Theorem 1 

Assuming GRH and LI, Rubinstein and Sarnak obtained an explicit formula for the 
Fourier transform of fiq-ai,...,ar in terms of the non-trivial zeros of Dirichlet L-functions 
attached to non-principal characters modulo g. More specifically they showed that 



(3.1) fig-au...,ar{ti,...,tr)=expii^Cg{aj)tj] Yl Yl Jo 



i=i / Xt^xo 7x>o \ + 7x 

X mod q 



First, we record an exponentially decreasing upper bound for fiq-ai,...,ar(^')y which is 
estabhshed in [6]. 
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Lemma 3.1 (Proposition 3.2 of [6j). Assume GRH and LI. Fix an integer r > 2. Let q 
be large and 0<e<l/2bea real number. Then, uniformly for all r -tuples of distinct 
reduced residues ( ) modulo q we have 

\i^q;ar,...,ar{ti, ■ ■ • ,ir)| < exp (-Cn (r )0(g) 1 1 1 1 1 ) , 

for ||t|| > 400 and 

\flq;au...,ar{tu • • • ,ir)| < exp(-Ci2 (r) £ ^(g) log 

for e < ||t|| < 400, where Cii(r) and Ci2{r) are positive numbers that depend only on r. 

Following the method developed by the author in [B] , we shall derive an asymptotic 
formula for fig-ai,...,ar (^ll^ll > AA/Var(g)^ in the range < A < -^/log log q, from which 
Theorem 1 will be deduced. 

Theorem 3.2. Assume GRH and LL Fix an integer r > 2. Letq be large andai, . . . ,ar 
be distinct reduced residue classes modulo q. Then in the range < A < A/log log q we 
have 

/i,;a„...,a. (||x|| > = {27:^/' [ CXp f + ' ■ ■ + 

^ ' -'l|x||>A V ^ / 

where 

F,4\) = (27r)-/2 / (x? - l)ixl - 1) exp f-^^^-^) d^. 

Proof of Theorem 1. First, it follows from Corollary 5.4 of [6J that 

max \Bg{a,b)\ < cf>{q). 

{a,b)eA{q) 

On the other hand. Proposition 5.1 of [6] yields 

\Bg{a, -a)| > 0(g), 
for all (a, g) = 1. Hence, we deduce 

(3.2) ma. ^ x 

(a,b)€A(q) Var(g) logg 

We also remark that this last estimate follows implicitly from the work of D. Fiorilli 
and G. Martin [2]. 

The first part of Theorem 1 follows from Theorem 3.2 upon using (3.2) and noting 
that Fj-fc(A) <^r 1. Let 1 < j < A; < r. If 1/4 < A < 3/4, then 

(2vr)-/2 / (a:? - l){xl - 1) exp f-^^^-^) dx > 5„ 
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for some positive number 5r which depends only on r. Moreover, we have 
(27r)-/^ / (X? - l){xl - 1) exp f-d±l_±A) ^ q, 

since the antiderivative of (x^ — l)e~^^/^ is —xe~^^^^. Hence we deduce, in the range 
1/4 < A < 3/4, that 

(3.3) ^j,fc(A) < -Sr, for all 1 < J < A; < r. 

On the other hand, it follows from (3.2) that there exist distinct reduced residues Oi, a2 
modulo q, such that 

(34) B,ia,,a2r J_ 

Var(g)2 ^"^log^g" 

Let as, . . . ,ar be distinct reduced residues modulo q that are different from ai and a2- 
Appeahng to Theorem 3.2 along with (3.3) we get 

Pi,;a,,...,a. (||x|| > XV^^) " i^^r^' [ exp ( -^i±—±^) dx 

S^oi^s)^ (loglogg)' 



2Var(5)2 (logg)3 

for some Kr > (which depends only on r), ii q is sufficiently large. Combining this 
inequality with (3.3) and (3.4) completes the proof. □ 

The first step in the proof of Theorem 3.2 consists in using the explicit formula (3.1) 
to approximate the Fourier transform fiq;ai,...,ar iti7 ■ ■ ■ i^r) by a multivariate Gaussian 
in the range ||t|| <^ A/log g/Var(g). To lighten the notation, we set jj, = lJiq-ai,...,ar 
throughout the remaining part of this section. 

Lemma 3.3. Assume GRH and LI. Fix an integer r > 2 and a real number A — 
A{r) > 0. //||t|| < Av/logg, then 

tr \ f tl + --- + t 



" I vmr ■ • ) ^ '''' i-^^j y-^) 

Bq{ajj^, aki)Bq{aj^, ak^)tjj,j^tk^tk2 + Or 



2Var(g)2 ^^.^^^^ .v^..,^...^..^..^..^.. V(log^)^ 

Proof. Using that Var(g) ~ 0(g) logg and |Cg(a)| — q°^^\ we infer from (3.1) that 



ELix(%)^j 



X mod q ' 
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Moreover, since log Jo{s) = -sV4+0(s^) for \s\ < 1, and ^ ^^^o T.y^>o l/(i+7x)^ < 

X mod q 



2Var(g), then 
(3.5) 



.2 



X mod g 



in the range ||t|| < A^/\ogq. Furthermore, the main term on the RHS of (3.5) equals 



(3.6) ^'""'^^ 



+ ■ ■ ■ + 1; 



2 Var(g) 

i<j<fc<r 

On the other hand, in our range of t, we have 



;73 Y Bg{aj,ak)tjtk. 



\ l<j<k<r / l<j<k<r 

which follows from (3.2). Combining this estimate with (3.5) and (3.6) completes the 
proof. □ 

Let $(a;) = e"^^/^ and denote by the ra-th derivative of $. Then $'^-'^^(2;) = 
— xe~^^/^, $*^^^(x) = (x^ — l)e~^'^/^, and more generally we know that $*^"^(x) = 
(— l)"iJ„(a;)e~^'^/^ where Hn is the n-th Hermite polynomial. We record the follow- 
ing result, which corresponds to Lemma 4.2 of [6]. 

Lemma 3.4. Let ni, . . . ,nr be fixed non-negative integers, and M be a large positive 
number. Then for any (xi, . . . , Xr) € W, we have 

Our last ingredient to the proof of Theorem 3.2 is the following basic lemma 

Lemma 3.5. If f\, . . . , fr are real valued functions, such that fi is odd for some 1 < 
i < r then 

[ fliXl) ■ ■ ■ friXr) exp ( - ^^ ~^ "^^'' ^rfx = 0, 

JVi<\\x\\<V2 \ ^ / 

for aUV2 >Vi>0. 
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Proof. Using the change of variables i/j = Xj ii j ^ i and Ui — —Xi, we deduce that the 
integral we seek to evaluate equals 

- / hivi) ■ ■ ■ UVr) exp f- ^'^V^^' ) dy, 

•>'Vi<||j/||<y2 V ^ / 

since /j is odd, which establishes the lemma. □ 



Proof of Theorem 3.2. Let R :— •\/Var(g) log logg. First, appeahng to Theorem 2 we 
get 

(3.7) 1^ (||y|| > AVVar(g)) = (AVVar(g) < ||y|| < R^+0 (exp (-(log log g)^/^^) , 

if q is sufficiently large. Next, we apply the Fourier inversion formula to the measure 
II, which yields 

/.(AVVa^<l|y|l<^) = (27r)-^ / / 6''-''^'+-+'^^^^ fi{s,, . . . , Sr)dsdy. 

^ ' JAi/Var(q)<||y||<R J s&T 

Let A = A{r) > r be a suitably large constant. Then, using Lemma 3.1 with e :— 
yl(Var(g))^^/^\/log q we derive 

f e'^''y'+-+'^y^^fi(si, Sr)ds = f e'(^ifi+-+^'-s"-)/i(si, . . . , Sr)ds + O ( ^) . 

Jsm^ J\\s\\<e \Q J 

Collecting the above estimates gives 

ix{\\y\\> A/VM^) =(27r)~'' / / e^(^i^^+-+^'-2"^)/i(si, . . . , s,)c^sdy 

+ 0(exp(-(loglogg)=^/2)), 
using that K^q^^^ <^ g^''. Upon making the change of variables 



tj :— •\/Var(g)sj, and Xj :- 



\/Var(Q) ' 

we infer that /i ^||x|| > A-\/Var(g)^ equals 
(3.8) 

(27r)-'' / / ^iit,,,+...+u.r)f^ ( ti ^ ^ tr \ ^^^^ 

A<||x||<iogioggi||t||<AVK^ \VVar(g) •,/Var(g) y 

+ 0(exp(-(loglogg)3/2)), 

Now we use the asymptotic expansion of fi (tiVar(g)~^/^, . . . , trVar(g)~^/^) proved in 
Lemma 3.3. First, the contribution of the error term in this asymptotic to the integral 
in (3.8) is 

^ (log log g)'" 
(logg)3 
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Next we shall compute the contribution of the main terms. Appealing to Lemma 3.4, 
we obtain 

(27r)-'^ / f g»(tixi+-+t.x.) exp f J^^ + " ' + ^r \ ^^^^ 

Jx<\\:>c\\<\oglogqJ\\t\\<AV^ V ^ ) 

- (27r)-'-/2 / exp (j'\±l_±^ + O f ^"i 

A<iix||<iogiog(? V 2 y V? / 

= (27r)-'-/^ / expf-^i±— ±^')dx + 0(exp(-(loglogg)3/^)), 

if q is large enough. Similarly, we deduce from Lemma 3.4 that for 1 < j < A; < r we 
have 

(27r)-'- / / i,ifce^(*i^i+-+*'^^'-) exp f- ^' ^ ^ dU^ 

A<||x||<loglogq J||t||<AVl5i^ \ 2 / 



dtd'x 



= -(27r)-''/2 / x.Xfeexp ( - ^^ ^ ^ )dy: + 0[^] =0 

A<||x||<loglogg 

which follows from Lemma 3.5. Furthermore, if 1 < ji < A;i < r and 1 < j2 < ^2 ^ f, 
then a similar argument along with Lemmas 3.4 and 3.5 shows that 

{2nr [ [ t,,t.,t,.t..e^(*-^+-+*-^) exp f- ^' + ^ + ^H 

^A<||x||<loglogq^||t||<Av15i^ V ^ / 

is < q-^ if ji 7^ j2 or ki ^ k2, and equals Fj>(A) + O (exp (-(log log g)^/^)) , if 
ji — 32 — j and ki = k2 = k. The theorem now follows upon collecting the above 
estimates and using Lemma 3.3. □ 

4. Very large deviations of iiq-ai,...,ar- Proof of Theorem 4 

In this section we prove a precise estimate for ^q-ai,...,ar{\x\oo > V) in the range 
V / {(l){q) log^ q) — 7- oo as g — 7> oo. The advantage of using the maximum norm is that 
the bounds in Lemma 2.1 can be made sharp, so that a precise estimate for large 
deviations of ^iq-ai,...,ar would follow from a close investigation of the tails of Y{q). 
Indeed, we have 

Lemma 4.1. Let e > be small. If q is sufficiently large, 2 < r < 0(g) — 1 and 
V > ■\/(t>{q), then for all distinct reduced residues Oi, . . . , modulo q, we have 

P (y(g) >V + q')< //,;ai,...,a.(k|oo > 1^) < 2rP (y(g) > ^ - g^) . 

Proof. The result can be derived along the same lines as in the proof of Lemma 2.1 
upon noting that 

r r 

f^r,au...Moo >V)< ^P(|X(g,a,)| > V) < {\Y{q)\ > V - \Cq{aj)\) , 

j=i j=i 

and \Cg(a)\ = g''^^) for all (a, g) = 1. □ 
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In view of this lemma, it suffices to obtain tfie analogous estimate for 

p,{V) ■.^V{Y{q)>V). 
Since E(exp(i cos(27r^(7^)))) = 7o(^), then 

(4.1) C{s) := r se'"" p,{V)dV = H 11 ( ' all s > 0. 

X mod q ' 

To prove Theorem 4, we shall estabhsh an asymptotic formula for log£(s) in a wide 
range of s, and then use the saddle-point method to extract the desired estimate for 
Pqiy) from that of His). We first collect some useful estimates and properties of the 
Bessel function /o(^)- 

Lemma 4.2. log/o(t) is a smooth function with bounded derivative on [0, +00) and 
satisfies 

'0(t2) ifO<t<l 
t + 0{log3t) ift>l. 



log/oW 



Proof. The first estimate follows from the Taylor series Io{t) — Y1^=q{^/'^Y^ 1'^'^-^ ■ More- 
over, we have 

(4.2) hit) = f e'^'^^^^'^de <e\ 

Jo 

On the other hand, taking e = 2^ we deduce 



Io{t) > [ e*'°'^^''^U9 > ee*^°^(2"^) > ^* 
Jo 



IOtt^ 



This together with (4.2) yields the second estimate. Finally, since Io{t) is a positive 
smooth function on [0, -|-oo) then log/o(t) is smooth and 

/g^cos(27r^)e*=°^(27re)^^ 



iaog/o(t))1 = 



< 1. 



□ 



Let f{t) be the real valued function defined by 
We prove 



logJo(t) ifO<t<l 
loglo{t)-t iit>l. 



Proposition 4.3. Iflogs/logq — >■ 00 as g —>■ 00, then 

£{s) =exp (^^^ -s log^ s + ^^slog s + 0{s(l){q) log q) 

V 27r TT 
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where 



\ p\<i j 

Proof. First, it follows from Lemma 4.2 that log/o is a Lipshitz function. Therefore, 
using (4.1) and Lemma 2.4 we get 

(4.3) iog£(.)= Yl Ei°g^o(^i^)= E El°g^o(^)+^2, 

X niod ' V mnH 



X mod 5 



where 



Xt^XO 0<7x<2W4+7x Xt^XO 7x>2 

X mod g " ^ mnH « 



X mod q 



<siV,(2) + s E E I^T 



Xt^XO 7x>2 4 
X mod q 



+ 7^ 



Furthermore, note that 

(4.4) y: E'°s'o(^) = 2. E E 7+ E E/(7)- 



Xt^XO 7x>2 
X mod q 



Xt^Xo 2<7x<2s 
X mod q 



Xt^Xo 7x>2 
X mod q 



Now, we infer from Lemma 2.4 that 



E E 



Xt^XO 2<7^<2s 
X mod g 

_ (0(g) - 1) 

27r 



2s 



+ 



^2 



2s 



logs + ^ f^*^^^) " + ^) + ^^'^^^^ 



iog2 5 + + (0(5) _ i)(i + log 2)) logs + O(0(g) logg). 



On the other hand, the second sum on the RHS of (4.4) equals 

r/(i)-.(*)^-r/(iy(^*>os. 



where 



s 

Es^\f{s)N,{2)\+<P{q) - 



\og{qt)dt -C s0(g) logg. 



which follows from Lemmas 2.4 and 4.2 along with the fact that \\m.t^^ f (y) Nq{t) 
0. This yields 



f /(f).^.(.) 



2s^ f<P{q)-i^^^^^m + <PiQ)-i 



27r 



27r 



dt+0{s(p{q) log q). 
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Making the change of variable u — 2s /t, we deduce that the integral on the RHS of 
the previous estimate equals 



^slogs r ^du + 0(s(j)( q) log q) 

- Mdu + 0{s4>{q)logq), 



TT Jo U 

using that R{q) ~ 0(g) logg by Lemma 2.3 and 

/(M)logM, ^ r f{u) , Plogw, logs 







which follow from Lemma 4.2. The proposition follows upon collecting the above esti- 
mates. □ 

Proof of Theorem 4- Assume that V is such that V/ {(p{q) log^ g) — >■ oo as g — >■ oo, and 
let s — s{V) > 1 be the unique solution to the equation 

(4.6) 2TrV = (0(g) - 1) log' s + 2 {D{q) + (t){q) - 1) log s. 

Then 



, D{q) ^ { D{q) V 27rV 
(4.7) g = exp - - 1 + t/ ( ,/ X -. + 1 + 



0(g) -1 VV0(9)-1 ; 0(g) -1 

Let e,S > he small real numbers to be chosen later and define 

si := (1 + e)s, S2 = (1 - e)s, and Vi := (1 + S)V, V2 := (1 - 5)V. 
Moreover, let 

poo pV2 

Ii := / se^^pq{t)dt, and h := / se"^ Pq{t)dt. 
Jvi J-00 

Then 

/"°° ,-"to ^si*^ ^ „-6gVi -^('^1) 



£(s) (1 + e)£(s) Jv-i (l + e)£(s) 

Therefore, appealing to Proposition 4.3 we obtain that the RHS of the previous in- 
equality equals 

exp ^-esV^i + ^^^2^ ^ (gilog'gi - slog's) + ^^^(silogsi - slogs) O(s0(g) logg) ) , 
which in view of (4.6) becomes 

exp {^-e5sV + ^^^|^slogs((l + e) log(l + e) - e) + O(s0(g) logg)^ . 
Combining this estimate with (4.8) and using that (1 + e) log(l -|-e) — e < e', we deduce 
<exp (-e5sV + e' ^^^^ ~ ^ s log s + Q(s0(g) log g) ) . 



£(s) ~ V 
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Now, (4.6) yields 2ttV > (0(g) — 1) log^ s. Therefore, choosing 6 = At/ log s wc get 



C{s) 



- (^l^*^*^^^^ ~ I)slogs + O(s0(?)logg)^ 



Furthermore, taking e = C-^log q/ logs for a suitably large constant C > 0, we obtain 

h 



A similar argument shows that 

h 



< exp(-s0(g) log?) . 



< exp(-s0(5) logq) . 



I 



Combining these two inequalities with Proposition 4.3 we deduce 

rV(i+s) /(h(a) - 1 D(a) \ 

/ se'^pq{t)dt = exp slog^sH — slogs + 0{s(j){q) logq) . 

Jv{i-S) V 27r TT J 

On the other hand, since V <^ (l){q) log^ s then 

se^'^dt = exp IsV + Oi s4>{q) V log s log g ) ) . 

'V(l-5) ^ ^ 

Hence, using that pq is a non-increasing function we infer from the two previous esti- 
mates that 

(4.9) p^{V{\^S)) < exp (^- '^^^^^~ ^ slogs + O (s0(g) ^logglogg)^ < Pg{V{l-S)). 
Now, using that D{q) ~ (f){q) logg by Lemma 2.3 we get 

and 



logg / 0(g)log^g^ 



y logs V ^ 

Moreover, we have 
(4.12) 



exp 



0(g) - 1 V v0(?) - 1 y - 1 

Finally, the theorem follows upon using (4.6) and inserting the estimates (4.10)-(4.12) 
into (4.9). 

□ 



20 
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